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Abstract 

We study the robustness of geometric phase in the presence of parametric noise. For that purpose 
we consider a simple case study, namely a semiclassical particle which moves adiabatically along a 
closed loop in a static magnetic field acquiring the Dirac phase. Parametric noise comes from the 
interaction with a classical environment which adds a Brownian component to the path followed by 
the particle. After defining a gauge invariant Dirac phase, we discuss the first and second moments of 
the distribution of the Dirac phase angle coming from the noisy trajectory. 

PACS: 03.65. Vf, 02.50.Ey, 03.67.-a 

1 Introduction 

The first reference to the role played by geometric phases in physics dates back to the work of S. Pan- 
charatnam [1] in the context of interferometry of polarized beams of light. Later the same phenomenon 
was described by M. V. Berry |f2l for quantum mechanical systems in the adiabatic limit. A mathe- 
matical insight into its origin was provided by B. Simon [31 which recognized that Berry phases could 
be interpreted as holonomies on a fiber bundle. Subsequently, quantum geometric phases (i.e. quantum 
holonomies) have been predicted and observed in various physical systems, and several generalizations 
and extensions were proposed flU 13 (see Q and the references therein). Generally speaking, we 
may say that geometric phases appear in correspondence with a cyclic evolution of a relevant Hilbert 
(sub)space. The dimension of the cyclic (sub)space determines the features of the corresponding geomet- 
ric phase: Abelian, i.e. U(l)-valued, for Hilbert space of unit dimension, nonAbelian, i.e. U(A r )-valued, 
in the case of A^-dimensional cyclic Hilbert space. 

Few years after the seminal papers by M. V. Berry and B. Simon, the scientific community discov- 
ered the potentialities of quantum mechanical systems in the context of information and communication 
technology JSJ, leading to the birth of quantum information science @- Since information processing 
obeys physical laws, a reversible quantum algorithm is represented by a unitary transformation as long as 
information is encoded as vectors in a Hilbert space. It follows that quantum holonomies, being unitary 
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transformations, can serve as quantum logical gates to implement quantum algorithms. This idea was 
first proposed and discussed in iflOl . where the authors also proved that universal computation ifTTt can be 
in general realized by means of solely (nonAbelian) geometric phases. A vast literature followed, which 
included both theoretical proposals 11121 [T3l and experimental realizations lfl4l [131 of simple holonomic 
quantum gates. 

In view of the application for quantum information processing, one usually considers the case of 
adiabatic geometric phase (also called Berry phase). In this case the system Hamiltonian is supposed to 
be a smooth function of the coordinates on a suitable manifold, often called the 'parameter manifold'. 
The geometric phase arises in the adiabatic approximation in correspondence with a closed loop in the 
parameter manifold. 

It is worth remarking, however, that that geometric approach to the computation can be rather de- 
manding from a technological point of view. Nevertheless, its advantage with respect to the standard 
dynamical approaches relies on the fact that Berry phases are argued to be particular robust with respect 
to noise. In particular, it is argued that holonomic quantum gates can be robust with respect to certain 
kind of classical parametric noise ifToll . This kind of noise can be modeled as coming from the interaction 
with a classical environment. 

The robustness of quantum logic gates is indeed a crucial issue because of the inherent fragility of 
quantum mechanical systems. Hence, although quantum error correction protocols exist fTTl . quantum 
gates which are a priori robust are welcome. It is worthwhile to mention another remarkable proposal for 
a fault tolerant computation, namely topological computation lfl"8l . This approach is based on nonAbelian 
Aharonov-Bohm topological phases. 

As we have anticipated, the geometric phase is believed to be robust with respect to classical para- 
metric noise. This argument has been the central issue of several investigations, in particular we mention 
the work by the G. De Chiara and M. G. Palma |fl9l , as well as the results presented in |[20l . Several 
physical models have been taken into consideration to the study of the robustness of geometric phase. 
Recently, the robustness of geometric phase has been tested experimentally in trapped polarized ultra- 
cold neutrons [21]. However it is worth noticing that this issue is largely independent of the details of 
the physical model under consideration. For this reason, our discussion will be devoted to the simplest 
settings in which geometric phases appear. In the following sections we consider a charged semiclassical 
particle which is adiabatically and cyclically moved in a static magnetic field. The particle acquires the 
Dirac phase which is proportional to the magnetic flux enclosed by the particle trajectory. Two different 
settings will be described: in the first case the semiclassical particle is in the presence of a homoge- 
neous magnetic field; in the second example the particle is subject to the field generated by a magnetic 
monopole. 

2 Dirac phase under noise 

In the following sections we consider the effects of parametric noise on quantum holonomies in two 
simple but remarkable examples. We consider as well two models for the noise, respectively represented 
by a Wiener and by a Ornstein-Uhlenbeck process. The examples involve a simple physical system 
made of a charged semiclassical particle in the presence of a static magnetic field B = (B x , B y: B z ). We 
indicate the corresponding vector potential as A = (A x , A y ,A z ) and the particle position as r = (x, y, z). 

Let us initially consider the case in which the particle follows a closed loop during a certain opera- 
tional time T, hence we have tq(T) = ro(0). In the adiabatic approximation, if the internal state of the 
particle in initially described by the a state vector *(0), after the closed loop in the parameter space the 
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internal state of the particle will be described by the vector ty(T) = e l ^ s \I/(0). Where 

4>g = -\ j> r A 'dr (1) 

is the acquired Dirac phase. Here we neglect the dynamical contribution to the phase factor. Let now the 
particle be subjected to Brownian motion due to the interaction with a classical environment. In this case 
a stochastic component r n (i) adds to the drift motion leading to the noisy trajectory 

r(t) = r (t)+r n (t). (2) 

Our aim is to study the Dirac phase acquired as a consequence of the noisy trajectory. First of all, we 
need a well defined notion of Dirac phase for the Brownian trajectory. The point is that the noisy path 
is in general noncyclic, while gauge invariance of the phase £T|) requires a closed loop. As we discuss 
below, it is possible to define a gauge invariant Dirac phase in a rather natural way. Secondly, since the 
Dirac phase is determined by the stochastic trajectory of the particle, we expect that the resulting phase 
is a stochastic variable as well. Our interest will be focalized on its mean value and variance. 

To conclude this section, we introduce our definition of Dirac phase for noncyclic trajectories. From a 
general point of view, it is possible to define the geometric phase in correspondence of noncyclic evolution 
if one introduces a rule that allows to close an open loop. As it was discussed in j6j, the natural choice 
is to close the loop with a geodesic curve, where the metric is defined by the hermitian product in the 
Hilbert space which is relevant in the given context. In our case, the natural choice is to consider the 
Euclidean metric in R 3 , hence we consider the following definition of gauge invariant Dirac phase: 

<j) g = J A ■ dr + J A ■ dr, 

where G indicates the straight line joining the final point r(T) to the initial point r(0). 



3 Particle in homogenous magnetic field 

In this section we consider the case of a static homogeneous magnetic field B = (0, 0, B). We can write 
the vector potential in the asymmetric gauge as A = (—yB, 0, 0), hence the Dirac phase is determined 
by the integral 

(j> g = f y(t)dx(t). (3) 
J o 

We are going to consider a noisy path of the form with a trivial drift component 

r (i) = (x(0),y(0),z(0)) 

yielding a trivial noiseless Dirac phase. For the sake of simplicity we consider the motion of the semi- 
classical particle as confined in the plane x — y, perpendicular to the magnetic field. 



3.1 Wiener process 

As a first model for the noise component, we consider the case of a Wiener process. Hence we impose the 
following conditions on two-times correlation functions of the Brownian component r n (i) = (x(t), y(t)): 

(x(s)x(t)) = B x 6(s-t), (4) 
(y(s)y(t)) = B y 5(s-t), (5) 
(x(s)y(t)) = 0. (6) 
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The noisy path has initial point x(0), y(0) and final point x(T), y(T). In order to ensure gauge invariance, 
we add one term to the equation (O, obtaining the following expression for the gauge invariant Dirac 
phase angle: 

4)9 = So y{t)dx{t) + \ [X(0) " X(T)] [y(0) + y(T)] = I y ( t)dx ^ + \ AS ' (7) 

where A := x(0) — x(T) and £ := y(0) + y(T). With this definition we can compute the corresponding 
mean value and variance of the phase factor e 1 ^ 3 . However, for small perturbation we can consider 
directly the mean and variance of the angle 4> g . Concerning the mean value, it is immediate to see that it 
vanishes as long as the two processes x(t) and y(t) are statistically independent. On the other hand, the 
variance reads 

°l = (f y{t)dx{t) j\{s)dx(s)) + (EA £ y{t)dx{t)) + \{Z 2 A 2 ). 
From the relations dU)-©, we obtain the following expression: 

a 2 g = \B x B y T 2 . (8) 
Hence we obtain that the variance of the Dirac phase angle 

grows linearly with the operational time T. 



3.2 Ornstein-Uhlenbeck process 

In the discussion concerning the robustness of geometric phases, a crucial role is played by a typical time 
scale characterizing the noise. A time scale is as well needed in order to deal with the adiabatic limit. 
For these reasons we consider a model of colored noise which presents a typical time scale, namely the 
Ornstein-Uhlenbeck stochastic process. Within this model, the coordinates of noise component fulfill the 
following stochastic differential equations: 

dx{t) = -T x x{t)dt + ^/D~ x dW x {t), 
dy(t) = -T y y{t)dt + ^/TT y dW y {t), 

where W x (t) and W y (t) are two independent normalized Wiener processes satisfying 

(dW x (s)dW x (t)) = {dW y {s)dW y {t)) = 5{s - t), 

and 

(dW X ( S )dWy(t)} = 0. 

The two-times correlation function decays exponentially: 

(x(s)x(t)) = 4e- r -l-*l, 
(y(s)y(t)) = e^e-r.l-*!, 
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where e 2 = e 2 = ^ and (x(s)y(t)) = 0. 

The Dirac phase angle can be written as the mean square limit (see e.g. |[22l ) of the following quantity 

4> g = lim S + -£A, 

AJ — >oo Z 

where 

Af-i 

S = y(tj)(x(tj + St) - x(tj)), 
3=0 

with 5t = T/M, and the proper term has been added to ensure gauge invariance. 
The the variance of the Dirac phase angle is given by 



a 2 = lim(5 + -SA) 2 = lim(cS 2 ) + lim(SEA) + 7 (£ 2 )(A 2 ). (9) 



^SA) 2 = lim(S 2 ) + lim(cSSA) + - 
The first term on the right hand side of equation Q is the limit of 

52 = ED^W*;))^ + 6t ^ -x(ti))(x(tj + 5t) -x(tj))), 

i 3 

where the average of y and x factorizes for the statistical independence of the processes. We have 

AS 2 = ^2(y(ti)y(t,)) [(x(t i+1 )x(t j+1 )) - + {x{U)x{tj)) - (x^x^-))] 

ij 

Evaluating the two-times correlation functions, and putting D x = D y and F x = T y ,one obtains: 

AS 2 = e 2 ^(yfcMtj)) [2e- r ^ st - e - T ^ +1 \ 5t - e -r|i-<+H*t} . 

ij 

The term in curled brackets is 

2(1 - e~ rst ) ~ 2T5t for |t - j\ = 

-(1 - e~ m ) 2 ~ -F 2 {8t) 2 for \i - j\ = 1 , 

e -r|wl«(2 - e ™ _ e - r5t ) ~ e'^^r 2 ^) 2 for |i - j| > 1 

and (y(ti)y(tj)) = e 2 e~ r ^~ J ^ t . Taking the limit <5i — > 0, only the terms with \i — j\ = do not vanish, 
leading to 
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V 2T5t ~ e 4 / 2rdt = 2e 4 TT = 2e 4 iV, 
. Jo 



j 

where N := TT is interpreted as the average number of statistically independent fluctuations. 
The second term on the right hand side of (O is the limit of the quantity 

<(x(0) - x(T))(y(0) + y(T)) £ yfoXxfc + - x(t,))), 

j 

which equals 
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and, in the limit of 5t — > reads: 

-r e 4 ^(e- r ^ +e- rT e ri O(e- r ^ +e- rT e r ^)^ 



-4Te 4 e~ rT 



f T , 








/ cosh 




cosh 


KH] 


Jo 









dt. 



The last integral reads 



-2e 4 r 



1 - e 



-2PT 



2r 



+ T 



Finally, the third term on the right hand side of © is 

Summing all the contribution, and taking the limit e~ rT — > 0, we can write: 

a 2 ~ 2e 4 (IT + 1) . 

From the last equation we see that, in contrast to the case of the Wiener process, for the Orstein-Uhlenbech 
process the variance 

a ~ eV 2 (rr+ 1). 

grows with the square root of the operational time. 



4 Particle in the field of a magnetic monopole 

In this section we consider another simple physical example, in which the semiclassical particle is sub- 
jected to the field of a magnetic monopole. If the magnetic monopole is sitting at the origin of the 
reference frame, we can write the corresponding vector potential, for z/R ^ —1, as follows: 



1 — z/R - , , , 

A • dr = — 7T {—ydx + xdy) , 

x z + y z 



(10) 



where R 2 = x 2 + y 2 + z 2 , or simply in spherical coordinates, for i? / it, as: 

. j 1 , 
A • dr = — tan —dtp, 
R 2 

where cos d = zj ' R and tan ip = y/x. For sufficiently small amplitude of the noise and short operational 
time we can take a linearized version of (fTOl : 

A-dr~A L -dr = - [fo + f x {x - x ) + f y (y - y ) + f z (z - z )] dx 
+ bo + 9x{x - x ) +g y (y - y ) + g z (z - z )j dy. 

In this approximation, we can write the following expression for the gauge invariant Dirac phase angle: 



+ 



A L ■ dr 



1 



f - -(f x A x + f y A y + f z A z 



A, 



9o ~ ^{g x A x + g y A y + g z A z 



A 



Vi 



(11) 



where A x = x(0) - x(T), A y = y(0) - y{T), A z = z{0) - z(T). 
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4.1 Wiener process 



In this section we consider a tri-dimensional Wiener process as noise model. We consider a trivial noise- 
less loop to which the Wiener process is superimposed. Denoting as B x , B y and B z the diffusion con- 
stants, we obtain the following expression for the mean value of dTTT) : 

{<t> 9 ) = ^(Al) - |<A*) = l -(g y B y - f x B x )T, 

which does not vanish in general and grows linearly with the operational time. 
Regarding the corresponding variance, we obtain: 

°i = m - {<t> 9 ? = \bt\ 

Hence the variance 

grows linearly with the operational time, and 

B = B x (3B x f 2 + 2B y f 2 + 2B z f 2 z ) + B y (2B x g 2 x + 3B y g 2 y + 2B z g 2 z ). 

4.2 Ornstein-Uhlenbeck process 

In the examples discussed above we computed the variance of the Dirac phase angle caused by a Brownian 
motion of the semiclassical particle. We explicitly considered the case of a trivial noiseless loop, in which 
the trajectory is purely Brownian. We obtained that the leading term in the variance is of the second order 
in the amplitude of the noise. Moreover, it increases linearly with the operational time for the Wiener 
process, and with the square root in the case of the Ornstein-Uhlenbeck process. This behavior can be 
compared with the results presented in fT9l , where the variance of the Berry phase was computed in the 
presence of a noise component modeled by an Ornstein-Uhlenbeck process. This noisy component is 
superimposed to a drift loop which is a precession around the z axis. In that case it was shown that the 
leading term in the variance is of the first order in the amplitude of the noise. Moreover, it decreases 
linearly with the operational time, leading to negligible fluctuations in the Berry phase. 

It is interesting to compare the variances of the Dirac phase obtained in the case of different drift 
loops. We have considered both the case of a trivial noiseless component 

r (t) = (sin(?9o)cos(i^o),sin(?9o)sin(^o),cos (i? )), 

and the case of a precession about the z axis described by the loop 

ro(t) = (sin ($o) cos (ipo + 2nt), sin (i)q) sin (ipo + 2nt), cos ($o))- (12) 

We have numerically simulated (following ll23l ) an Ornstein-Uhlenbeck process affecting these loops 
and estimated the variance of the corresponding Dirac phase angle. The results are plotted in figure [T]for 

cos $0 = l/v 7 ^- 

We notice two different pattern of the variance of the Dirac phase angle as function of the average 
number of fluctuations in the noisy component. In the case of trivial drift loop (purely Brownian motion) 
the variance always increases as the square root of the average number of fluctuations. On the other hand, 
for nontrivial noiseless loop (Brownian component superimposed to precession) a transient behavior is 
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present in which the variance decreases with the number of fluctuations. This behavior is in agreement to 
what was found in fT9ll and is due to contribution of the first order in the noise amplitude. By increasing 
the value of N, the first order contributions become negligible while the second order ones become 
predominant. 




100 200 300 400 500 600 700 800 
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Figure 1 : The plot shows the numerically estimated variance of the geometric phase angle for a particle 
in a monopole field. The variance is plotted as a function of the average number of noise fluctuation 
for a Ornstein-Uhlenbeck process modeling the noise component, with amplitude e = 0.05. The data 
represented by dots refer to a trivial drift loop (purely Brownian motion), they are well fitted by a square 
root law o g = ay/N + b, with a ~ 0.0025 and b = —0.00016. The data represented by circles correspond 
to the drift loop in equation (fT2~l) (Brownian motion superimposed to precession). 



5 Conclusion 

We have computed the mean value and the variance of the Dirac phase acquired by a semiclassical particle 
subjected to Brownian motion. If the trajectory is purely Brownian the variance of the Dirac phase angle 
always increases as function of the operational time (or the average number of noise fluctuations). On the 
other hand, a transient behavior is observed if the Brownian motion is superimposed to a noiseless drift 
loop. 

In the case of pure Brownian motion, we have obtained an expression for the variance which is of the 
second order in the amplitude of noise and increases with the operational time. In particular, if the noise 
is modeled by a Ornstein-Uhlenbeck process the variance grows with the square root of the operational 
time. 

The case of the Dirac phase can be viewed as instance of geometric phase. Hence we can compare 
our results to others which refer to Berry phase. In lfl9l . it was shown that in the case of an adiabatic 
precession of a 1/2-spin the leading term in the variance of the Berry phase is of the first order in the 
amplitude of the noise. Moreover, this terms decrease linearly with the operational time. This behavior 
is in accordance with the transient behavior of the Dirac phase for nontrivial noiseless drift loop. Notice 
that the second order effects become relevant for long enough operational time. 
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At the best of our knowledge, the presentation of the effects of second order in the variance of the 
Dirac phase introduces a new element in the study of the robustness of geometric phases. We argue 
that second order effects are feasible to be observed in experimental settings as in l2TTl . Recently, the 
effects of non-adiabaticity in the noise component were studied in ll24l . The pattern of the variance of the 
corresponding geometric phase as function of the operational time is qualitative analogous to the results 
presented here, in the sense that the squared variance grows linearly in time. Quantitatively, this effect 
is of the first order in the noise amplitude in ll24l while in our analysis, which assumes the adiabatic 
approximation, the effect of the second order. 
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